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In recent experiments on thin 4 He films absorbed to rough surfaces Luhman and Hallock (Ref. 
attempted to observe KT features of the superfluid-normal transition of this strongly disordered 
2D bosonic system. It came as a surprise that while peak of dissipation was measured for a wide 
range of surface roughness there were no indications of the theoretically expected universal jump of 
the areal superfluid density for the strongly disordered samples. We test the hypothesis that this 
unusual behavior is a manifestation of finite-size effects by numerical study of the corresponding 2D 
bosonic model with strong diagonal disorder. We demonstrate that the discontinuous features of 
the underlying KT transition are severely smoothed out for finite system sizes (or finite frequency 
measurements) . We resolve the universal discontinuity of the areal superfluid density by fitting our 
data to the KT renormalization group equations for finite systems. In analogy to our simulations, we 
suggest that in experiments on strongly disordered 2D bosonic systems the very existence of the KT 
scenario can and should be revealed only from a proper finite-size scaling of the data (for 4 He films 
finite-size scaling can be effectively controlled by the scaling of finite frequency of measurements). 
We also show relevance of our conclusions for a wider class of systems, such as superconducting 
granular films, Josephson junction arrays, and ultracold atomic gases, where similar difficulties 
appear in experiments designed to verify KT transition (especially in disordered cases). 

PACS numbers: 68.15.+e, 67.70.+n, 64.70.Ja, 68.35.Ct, 67.40.Pm 



I. INTRODUCTION 

For decades 4 He films played a prominent role in 
the research of superfluidity of two-dimensional (2D) 
bosonic systems. The Kosterlitz-Thouless (KT) theory- 
describes the finite temperature superfluid phase transi- 
tion for the 4 He films absorbed to the flat substrates^ 
One of its main predictions is the universal jump from 
zero of the areal superfluid density ctsf at the transition 
temperature T cr given by the relation 4 
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where m is the mass of a 4 He atom. Theoretically, ac- 
cording to Harris criterion, disorder is irrelevant at the 
KT transition, and should not alter the discontinuous 
KT scenario^ Recently, Luhman and Hallock (Ref. [l]) 
performed extensive quartz crystal microbalance (QCM) 
experiments of the superfluid transition of the 4 He films 
absorbed to rough surfaces. Surprisingly, their obser- 
vations were not directly agreeing with the expected KT 
scenario of the 2D superfluid transition. They found that 
it was increasingly difficult to identify signatures of the 
universal jump JT]) as the surface roughness was increased 
(see also Ref. M). For strong enough disorder the pres- 
ence of the KT transition was seen only indirectly by the 
observation of a peak in the dissipation. Could their ob- 
servations be an indication of some strong-randomness 
scenario of the phase transition, or is there a way to rec- 
oncile those experiments with the standard KT theory? 

Note that for ideal (disorder free) 2D bosonic systems 
discontinuity of the superfluid density at the critical point 
is a characteristic of the phase transition in the thermo- 
dynamical limit of infinite system sizes. In the experi- 



ments with the superfluid helium films the sample size 
can simply set finite length scales, e.g., it can be the 
size of a confining channel^ or the radius of spherical 
grains R in the case of helium films absorbed in packed 
powders, 8 or the pore radius R in the case of helium 
films absorbed on porous glass^ (in the latter two cases 
crossover from the effective 2D to 3D behavior should 
be kept in mincU^). Finite length scales can also be 
related to finite frequency of measurements, when one 
is bounded by the frequency-dependent vortex diffusion 
length r w = y^lAD/oj, where D is the vortex diffusivity 
and lu is the experimental frequency . 11 ' 12 Finite-size ef- 
fects in experiments smear the universal jump ([T]), and 
disorder is expected to increase this finite-size broadening 
of the KT transition even more. We think this scenario is 
applicable to the finite frequency measurements on He 
films absorbed to the rough surfaces in Ref. The pur- 
pose of this study is to show that in these experiments 
disorder of the substrate can produce such dramatic ad- 
ditional broadening of the superfluid transition on top of 
finite-size effects so that it could be indeed difficult to re- 
solve the growth of the superfluid fraction on the varying 
normal fluid background. 

In Sec. |H] we start with the description of 4 He films 
absorbed to rough surfaces studied in Ref. [I| by a 2D 
bosonic model with diagonal disorder. We want to be as 
close in our simulations to the experimental situation of 
Ref. Ill as possible in order to guarantee that all physics 
that might be important for explanation of surprising 
outcome of Ref. [lj will be preserved. We represent the 
microscopically disordered rough surface by a random 
potential (i.e., diagonal disorder). We keep temperature 
fixed, and tune the superfluid transition by adjusting the 
chemical potential (the average helium film thickness is a 
monotonic increasing function of the chemical potential) . 



2 



In the ideal case, we observe the anticipated finite-size 
broadening of the universal jump of the superfluid density 
([1]) near the critical point of the KT transition. Inclu- 
sion of the strong disorder leads to such severe additional 
broadening that, compared to the ideal case, it may seem 
impossible to make any definite conclusions about the 
underlying transition type. Our analysis utilizes Monte 
Carlo simulations of different system sizes, or, equiva- 
lently, ability to set by carrying out measurements at 
different finite frequencies. We resolve the asymptotic 
universal jump of the superfluid density |T]) by showing 
that data do obey the KT renormalization group (RG) 
equation. Finally, once we have established the KT sce- 
nario in the strongly disordered case, in Sec. lIIII we carry 
out a qualitative comparison of our simulations with the 
aforementioned Ref. |l| QCM experiments with 4 He films 
absorbed to rough surfaces. In the context of the QCM 
technique, we also discuss how the dissipation peak is 
affected by disorder. 

With respect to the KT transition, we expect the com- 
bined effect of diagonal disorder and finite-size broaden- 
ing to be qualitatively similar for different physical sys- 
tems, and to be independent of the particular type of 
disorder realization and choice of a control parameter. In 
Sec. II VI we discuss how our results for the 4 He films can 
be related to superconducting granular films, Josephson 
junction arrays, and ultracold atomic gases. 



II. SIMULATIONS 

In Fig. [T] we show a phase diagram in the temperature- 
density plane for a 2D bosonic model ([2]) with strong 
diagonal disorder. The transition from the normal to 
superfluid phase happens according to the KT scenario 
at the phase boundary. At some finite temperature T we 
drive the system through the transition by increasing the 
boson density. At the critical point, the superfluid areal 
density experiences the universal jump <7sF C r according 
to p. 

For a system of a given finite size, the superfluid den- 
sity is expected to be a monotonic function of the boson 
density, and we study its sharpness with respect to disor- 
der. The strength of disorder also determines the position 
of the phase boundary. 

Note that the zero temperature quantum phase tran- 
sition is in a different universality class In the case of 
4 He films on disordered substrates, the threshold boson 
density uqm < 1 per lattice site corresponds to up to 
several atomic layers^ 



A. Model 

We start with the description of bosonic systems with 
diagonal disorder by a general lattice bosonic Hubbard 




Boson density 



FIG. 1: Generic phase diagram in the density-temperature 
plane for 2D bosonic system with strong diagonal disorder 
(?iqm is the threshold boson density). We study the superfluid 
transition along the dotted line path. 



model: 
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where a\ and hi are the creation and annihilation oper- 
ators, hi — a\a>i is the number operator on site i, and 
denotes summation over pairs of nearest neighbors, 
each pair counted once. The first term describes nearest 
neighbor hopping, and the second term corresponds to 
soft-core on-site repulsion. The chemical potential fi de- 
fines the average number of bosons per site n(ji), whereas 
Vi represents random on-site potential. The value of Vi 
is distributed uniformly between — A and A. 

We do not carry out quantum Monte Carlo simula- 
tions of the Hubbard model |(3J) in continuous time^. 
Instead, we map 2D quantum bosonic model ([2]) on the 
equivalent (2+l)D classical J-current modelJ^ For this 
classical (2+l)D counterpart of quantum 2D model (J2J) 
we perform classical Monte Carlo simulations with high- 
performance worm algorithmj 17 ' 18 Action of the so-called 
J-current model is written in terms of integer link-current 
variables J = (J x , J v , J r ) defined on a N x N x N T space 
and imaginary time lattice (N t St = h/k^,T, where St is 
an imaginary time interval): 



(^) 2 ]+|E[( J n) 2 -(A + W 



n n 
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where n = (r, r) are discrete space and imaginary time 
coordinates. J has to be conserved and therefore sat- 
isfies the local zero-divergence constraint, V • J n = 0. 
Spatial currents J x and J v represent hopping, and imag- 
inary time current component J T represents occupation 
number. We choose to consider only positive occupation 
numbers, J T > 0. This is a natural choice for atomic sys- 
tems, and it cannot affect the universal properties of the 
transition* 1 ^ Correspondence of parameters of the models 
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(J2j) and ([3]) are given by 
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where no is some offset boson density of model ©. 19 Here 
the value of v r is uniformly distributed between — A and 
A, where A = A/U. 

The total boson density a (proportional to the boson 
number per site n) and the superfluid density ctsf are 
determined from Monte Carlo estimators of the winding 
numbers Wi = l/Ni ^ n J^, and are given b y 16 i 20 
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where fh = ti 2 /2ta 2 is an effective boson mass on a lat- 
tice with spacing a. In 2D, the statistics of spatial wind- 
ing numbers are essentially discrete, thus an appropriate 
way to find K is through fitting collected winding num- 
ber histograms with a discrete gaussian P(W x ,W y ) oc 
cM-(W 2 + W 2 )/2K). 

For dcfinitcness, wc choose numerical values of param- 
eters to be close to those characteristic for experiments 
with He films at temperatures of about 1 K, i.e., hop- 
ping, interaction and temperature are considered to be 
of the same order t ~ U ~ k^T, and tlq ~ «qm ~ 1 to 
correlate with a transition of several layers of a 4 He film. 
Corresponding effective parameters of the (2+l)D clas- 
sical model ([3]) are t, and U. Our choice of finite L T is a 
compromise between the efficiency of discrete imaginary 
time codes and staying very close to the original model. 
We have varied N T , while t and U can be estimated from 
(j4|), so that the boson number density in the vicinity of 
the superfluid transition in the ideal system would stay 
about 1 per lattice site. We found t = 4, U = 0.2 and 
N T = 6 to be appropriate values. 



B. Results 

We model systems of finite spatial size N x N to be 
24x24, 48x48, 96x96 and 192 x 192. Three situations are 
considered: no disorder A = 0, and disorder of strengths 
A = 3 and A = 6. In the later cases, we average data 
over about 500 realizations of disorder. Temperature is 
kept fixed throughout all simulations, and we study how 
the average boson number per site n and the superfluid 
density ctsf behaves as we pass over the transition point 
increasing chemical potential jl, see Fig. [2] and [3] Note 
that it is more informative to study the dimensionlcss 
reduced superfluid density K rather than ctsf itself. Ac- 
cording to (|5b[) . the coefficient of proportionality relating 
K and ctsf is h 2 /ifi 2 k B T . Thus, the universal discontinu- 
ity of the superfluid density ([I]) at the fixed temperature 
KT transition is equivalent to the jump of the reduced 
superfluid density K from to K CI — 2/ir. 




FIG. 2: Average number of bosons per lattice site n as a func- 
tion of chemical potential jl. Three sets of curves represent 
the cases of no disorder A = 0, and disorder strengths A = 3 
and A = 6. Statistical errors are smaller than the symbol size. 
Data in each set for different system sizes 24 x 24, 48 x 48, 
96 x 96 and 192 x 192 collapse on one and the same curve, 
i.e., n is independent of system size. 



For a given disorder strength the absence of saturation 
of the reduced superfluid density K (/2, N) for different 
system sizes below the critical value K cr and collapse of 
the data above K CI , see Fig. [31 is a qualitative indication 
of the KT discontinuity at infinite system sizes. To per- 
form quantitative analysis we use the integral form of the 
KT RG equations (|8al8b|) (see Ref. [SI) 



dt 



K(ji,N 2 )/K CI 
K(Ji,Ni)/K cr t 2 (\nt-£,(fi))+t 



Here is the size-independent microscopic param- 

eter characterizing the vortex fugacity, and, in con- 
trast to the thermodynamical reduced superfluid den- 
sity K(fl) = K(Ji,N — > oo), £(/i) is an analytic function 
of jl (i.e., may be expanded into the Taylor series 
£(/}) w 1 + £'(0) • (ft. — fta) + ■ ■ ■ in the vicinity of the 
critical point /x cr ). Those two properties of £ we use to 
check the consistency of our data with the KT scenario 
of the transition. For different pairs of K(ft, N) we solve 
^ for £(//) and observe that its values collapse on the 
line near ft cr , see Fig. 01 Equation ^ indirectly repre- 
sents KT RG flow of K (ft, N) for increasing length scales 
Ni < N 2 < N 3 < N 4 < . . .. In the limit N — > oo it gives 
equation 



K(ft) 



(7) 



for the thermodynamical value of K(ft). Finally, using 
determined function £(/2), we solve ([7]) for K(fi), see 
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FIG. 3: Reduced superfluid density K — h 2 osf /m 2 ksT as a 
function of chemical potential fi. Three sets of curves repre- 
sent the case of no disorder A = 0, and disorder strengths 
A = 3 and A = 6. In each set, four curves from the bottom 
to the top correspond to finite system sizes 192 x 192, 96 x 96, 
48 x 48 and 24 x 24. Statistical errors are smaller than the 
symbol size, (a) The upper graph shows K in the wide range 
of p,, and (b) the lower graph depicts proximity to the KT 
transition, where the solid curves are the KT extrapolation 
to the infinite system size. Horizontal lines are equal to the 
value of the universal jump K CI — 2/n at the KT transition. 
The dotted curve on the upper graph (a) is an example of the 
extrapolation of the data to lager system size 10 4 x 10 4 . 



Fig. [HJd (solid curves). Note that, due to the logarith- 
mic form of the finite-size corrections to K near the KT 
transition^ K{N) — K cr oc 1/ InJV for p, > p, CI , data for 
systems of exponentially different sizes is required for the 
KT RG analysis. 

The critical observation, see Fig. [3^,, is that disorder 
dramatically increases finite-size smearing of the KT dis- 



FIG. 4: RG parameter £ as a function of p,. Three sets rep- 
resent the cases of no disorder A = 0, and disorder strengths 
A = 3 and A = 6. In each set, values of £ are deduced by RG 
analysis © of different pairs of K corresponding to pairs of 
systems 192 x 192 and 48 x 48, 96 x 96 and 24 x 24, 192 x 192 
and 24 x 24 at Fig. [3b. The solid lines are linear fits for all 
data points at a given disorder strength. 



continuity of the superfluid density ctsf oc K . For exam- 
ple, consider system of size 24 x 24. Compared to the 
ideal case (disorder strength A = 0), there is no clearly 
visible characteristic bending in the data trend for the re- 
duced superfluid density K at a strong disorder (A = 6) 
near the anticipated discontinuous KT transition. And, 
yet, the self-consistent RG analysis of the data for expo- 
nentially different system sizes will and is revealing the 
KT transition. 



III. COMPARISON WITH EXPERIMENTS 

Our simulations are directly applicable to explain the 
role of disorder in the superfluid transition in 4 He film 
absorbed to microscopically rough surfaces, see Ref. [l|. 
In these experiments, sample surfaces are prepared by 
coating QCMs with CaF 2 . The height distribution of a 
CaF 2 film is close to a Gaussian function, and its disper- 
sion increases with the nominal deposition thickness (the 
system of zero disorder is an uncoated plain QCM crystal 
with typical gold electrodes). The crystals are mounted 
in a cell. The QCM crystals are operated in the third 
harmonic /o « 15 MHz of their fundamental resonance 
frequency. At the constant temperature 1.672 K, 4 He 
is incrementally added to the cell, correspondingly, the 
thickness of the 4 He film absorbed to the sample surface 
grows. The resonance frequency of oscillations f = lu/2it 
of the QCM crystal with an absorbed 4 He film is directly 
related to the viscously clamped normal fraction of the 
4 He film. When the film undergoes a transition from the 
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normal to superfluid state, one expects to see a sharp 
step in the resonance frequency / because of the decou- 
pling of the zero-viscosity superfluid fraction of the 4 He 
film from the substrate^ The absorbed film shifts the 
resonance frequency from the value /o by —A/ which 
is proportional to the difference of total and superfluid 
areal densities of 4 He, —A/ oc a — osf ( see insert of 
Fig. [5]). At a finite (angular) frequency uj, one probes 
length scales of order of a diffusion length — yJ\Aj5juj, 
where uj ~ ujq — 27r/o, and D is diffusivity constant. 11 
This sets finite-size scaling for asp and smooths out the 
zero frequency sharp step in A/. It was found in Ref. 
that with higher disorder no indication of this step is ob- 
servable. Another signature of the superfluid transition 
in 4 He films is a peak in the dissipation, and it indirectly 
manifested the presence of the KT transition in these 
experiments. 

A. Important length scales 

The superfluid transition in 2D bosonic systems is 
driven by unbinding of thermally activated vortex- 
antivortex pairs. Kosterlitz and Thouless used the 
renormalization-group technique to consider the atten- 
uating effect of smaller pairs on the interaction between 
the respective members of larger pairs. They derived 
recursion relations for a scale-dependent superfluid den- 
sity osf (7) (or, corresponding reduced superfluid den- 
sity K(l) = itsf(0^ 2 /m 2 kBT) and vortex-pair excitation 
probability y 2 {l) 

dK- l {l)/dl = ATT 3 y 2 (l)+0[y i (l)}, (8a) 
dy 2 (l)/dl = 2[2 -TrK(l)]y 2 (I) + 0[y 4 (1)1 (8b) 

where / = ln(r/ao), and contribution of vortex- 
antivortex pairs of sizes less than r has been already 
taken into account. The vortex core size ao is of the or- 
der of the healing length of the superfluid, i.e., the length 
over which the superfluid density can change from zero 
to its mean value. Here, we explicitly note that KT RG 
equations are only valid for y 2 <C 1 (for y 2 3> 1 the vor- 
tex density is so high that it is ambiguous to group them 
into some hierarchy of vortex-antivortex pairs of different 
sizes and then to carry out renormalization). 

For an idealized experiment carried out at zero fre- 
quency on infinite size system, the recursion relations are 
iterated from ao to infinite scale. In the normal phase 
y 2 diverges at infinite scales, and <tsf is renormolized to 
zero. In contrast, the superfluid phase is characterized 
by vanishing y 2 and some finite ctsf- Phase transition 
from the superfluid to the normal phase corresponds to 
unbinding of vortex-antivortex pairs as they get excited 
to infinitely large scales of their mutual separation, and, 
in addition, the superfluid transition is accompanied by 
the universal jump |T]) of <7SF C r- 

In finite systems the maximum vortex separation is 
bounded by the system size, and the KT iterative renor- 
malization process is cut off at a corresponding scale. 



This broadens the sharp KT transition and leads to a 
finite-size rounding of the KT jump in ctsf, see Fig. [3] In- 
terpretation of experiments at finite frequencies requires 
the theory of Ref. [ll| that incorporates the dynamic re- 
sponse of the vortex plasma to an externally applied oscil- 
lating field. It was shown that vortex pairs of size larger 
than the vortex diffusion length r u = y/lAD/u> cannot 
equilibrate to the oscillating field of the frequency uj, and 
do not contribute to renormalization, so the iterations 
are truncated at r w . Thus, results of our simulations for 
asp for finite-size systems are in one-to-one correspon- 
dence to finite- frequency experiments, provided that uj is 
adjusted so that r u is equal to the simulated system size. 



Let's specify particular values of length scales for 4 He 
films absorbed to rough CaF 2 surfaces studied in Ref. 1. 
Average helium film thickness in the vicinity of the super- 
fluid transition in those experiments was reported to be 
about 4 layers ~ 15 A on the flat substrate (on the disor- 
dered CaF 2 substrate the amount of helium absorbed was 
larger). Temperature was in the range 1 — 2 K, thus the 
vortex core size is expected to be on the order of atomic 
dimensions, ao ~ 10 A. Because ao is comparable to the 
helium film thickness no 3D related vortex effects are as- 
sumed to exist. The diffusivity constant D was not mea- 
sured and might be affected by disorder, but from dimen- 
sional analysis D w h/m ~ 2 x 10~ 4 cm 2 /secii 1 - (A cau- 
tionary remark: for different substrates the numerical co- 
efficient before h/m in the formula for D was reported 1 - 
to be in the wide range from 0.1 to 20.) Dynamical mea- 
surements of (jgp were carried out in the course of op- 
eration of the QCM crystals at uj rs 2tt x 15 MHz, and 
this gives r w ~ 500 A. QCMs were of centimeter dimen- 
sions, thus their width is an irrelevant length scale. CaF2 
surfaces coating the QCMs consisted of peak-type struc- 
tures with an average separation R s < 100 A. Because 
the roughness scale is less than the diffusivity length, 
Rs *C r u , disorder is considered to be microscopic. In 
Ref. [J by varying deposition time, surfaces with aver- 
age height of the CaF2 structures from to ~ 50 A were 
prepared and used (0 stands for uncoated (plain) QCM 
crystals) . 



With respect to the lattice model ^ used in our sim- 
ulations in Sec. |TTJ the vortex core size ao stands for the 
lattice spacing a, and the ratio r u /a Q ~ 50 defines lattice 
spatial size to be N x N ~ 50 x 50. We consider height 
variations of the CaF2 profile to act as a random po- 
tential. We chose on-site potential disorder to represent 
the closely spaced, raised structures of CaF2. Particular 
choice of the disorder representation cannot affect univer- 
sal properties of the second order KT transition. How- 
ever, to be in line with experiments of Ref.[l| we must sat- 
isfy criterion of the microscopic disorder r w 3> R s - And 
indeed this is true for our lattice model, because it corre- 
sponds to R s — a ~ ao and system size r w /ao ~ N 3> 1 . 
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B. Broadening of the KT transitions 



The experiments of Ref. Q] with 4 He films absorbed to 
rough CaF 2 surfaces can be compared to our model in 
the following way (see (|5al5bp ^) 



— A/ cx a — ctsf an - jK, 



fha 2 



k B T/H 2 . (9) 



For our choice of parameters 7 ~ 1 (we choose 7 = 0.75). 
Our simulations, specifically Fig. [5j qualitatively resem- 
ble the experimental situation. Actually, Fig. [5] contains 
excessive information relative to those experiments of 
Ref. Q]. There, QCM measurements were carried out at a 
fixed frequency, i.e., r u was fixed, while CaF 2 substrates 
of different roughness were used. In other words, at a 
given disorder strength (i.e., on a given coated QCM) 
only a single curve for mass decoupling ^ as a function 
of the chemical potential was observed, for example, let's 
say that these were upper curves in each set on Fig. [5] 
for which r w corresponded to 192 x 192 lattice sites in 
our simulations. In full agreement with experimental ob- 
servations, we see that, as disorder strength is increased, 
finite-size behavior completely masks the step in A/ re- 
lated to the KT universal jump of the superfluid density 
osp. This is exactly the observation that did not give firm 
evidence to unambiguously claim the KT mechanism of 
the superfluid transition in Ref. Q]. If, at each disorder 
strength, QCM measurements were carried out at a se- 
ries of frequencies, then similar to our results, Fig. [5j 
one could have had enough data to perform the KT RG 
analysis. For further discussion see Sec. fVl 

An essential condition of the KT transition is an ob- 
servation of the dissipation peak during dynamic finite 
frequency measurements. It signifies change in the dy- 
namical response from the one determined by the free 
vortices moving diffusively at the normal side to the one 
produced by the vortex pairs of size greater r w failing to 
equilibrate to external oscillating field at the superfluid 
side. In Ref. Q] even for the strongly disordered samples 
the peak in dissipation was present and was the basis to 
claim the KT scenario of the superfluid transition. Our 
simulations are statistical mechanical and do not simu- 
late per se the dissipation peak of dynamical origin. How- 
ever, using the results of the linear dynamical theory for 
real and imaginary parts of effective dielectric constant 
of vortex plasma^ we can predict the shape of the dissi- 
pation peak from our numerical data. The dissipation is 
measured in terms of the change in quality factor Q that 
is given b y 12 ' 23 



oc 



(Ky) 2 



1 + {TT i Ky 2 f 



(10) 



The dissipation peak is positioned close to the point of 
the inflection of the growing ergp as one goes from the nor- 
mal to the superfluid phase while increasing the control 
parameter, chemical potential /i. But KT RG works in a 
narrow window of the transition, e.g., we had to use only 
a portion of the simulated superfluid density data Fig. \3jp 
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FIG. 5: Frequency shift theoretically defined —A/ oc n — 
yK as a function of chemical potential p,. The three sets of 
curves represent the cases of no disorder A = 0, and disorder 
strengths A = 3 and A = 6. In each set, the four curves 
from the bottom to the top corresponds to finite system sizes 
24 x 24, 48 x 48, 96 x 96 and 192 x 192. Statistical errors are 
smaller than the symbol size. The insert shows an idealized 
sketch of —Af(fi) in the vicinity of the superfluid transition 
at /i cr ; the step corresponds to the jump of superfluid density 
o"sf in the system of infinite size, see Q. 



to deduce from © the size-independent microscopic pa- 
rameter £. Too deep in the normal phase y 2 > 1, and 
the KT recursion relations (|8a|8b|) break down. Thus, 
we had to extrapolate our data to higher system sizes, 
so that y 2 < 1 and the dissipation peak is positioned in 
the same range of jl. For each disorder strength, we used 
linear fits for £(/!) from Fig. [4] to recalculate the finite- 
size smeared reduced superfluid density K{fi, N = 192) 
of the system of size 192 x 192 into K(ji, N = 10 4 )M The 
value of y 2 (jl,N) for known K(jx,N) and £(/i) is given 
by 



y 



1 

2^ 



k \k' 



(11) 



Our observations for A(l/Q), see Fig.[BJ are that for one 
and the same r w with an increase of disorder strength 
the dissipation peak will be broadened and will be shifted 
deeper in the normal phase away from the critical point 
of the superfluid transition. A high level of noise in the 
data for 4 He films coating rough CaF 2 surfaces, Ref. [l|, 
did not allow a validation of this statement (note that 
dissipation peaks were clearly seen for all values of dis- 
order). However, in those experiments there is room left 
to raise accuracy at least by the order of magnitude. 25 

We have also studied finite-size scaling of the width of 
the dissipation peak for a given disorder strength. Its 
broadening one-to-one resembles the finite-size broaden- 
ing of the superfluid density jump (correspondingly, po- 
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temperatures (in Ref. [l. T — 1.672 K) or, equivalently, at 
high average helium film thicknesses, see the phase di- 
agram Fig. [TJ For T > 1 K we expect that the surface 
tension predominates, and microscopic puddles of helium 
are created in the valleys of the disordered surface profile 
while the helium film thins close to the ridges of the peak- 
type structures of the CaF2 substrate. A study similar to 
Ref. [l] of the superfiuid transition of 4 He films absorbed 
to the same rough CaF 2 substrates was carried out at low 
temperatures T < 1 K. It turned out that in this case 
broadening of the KT transition was almost independent 
of the strength of disorder^ This allows us to speculate 
that a different regime from the one considered in this 
paper was achieved for T < 1 K, when van der Waals' 
forces win over the surface tension and the helium film is 
uniformly flattened over the rough surface, and, hence, 
the surface potential is uniform. 



FIG. 6: Dotted lines are the reduced superfiuid density 
K <x ersF and solid lines are the dissipation peak A(l/Q) 
as a function of chemical potential fi for the systems of size 
10 4 x 10 4 . The three pairs of curves represent the cases of no 
disorder A = 0, and disorder strengths A = 3 and A = 6. 
Chemical potential is shifted by the critical value p, CT of the 
KT transitions in the corresponding systems. 

sition of the peak shifts to higher chemical potential as 
r w is increased). 

C. Disorder description in detail 

Our simulations explain broadening of the KT tran- 
sition in microscopically disordered 2D bosonic systems 
such as the 4 He films absorbed to rough CaF 2 substrates 
m Ref. 1. In our case, R s <C r u , the roughness scale 
R s is much less compared with the dimension r u of the 
region effectively studied at finite frequency. Different 
case of macroscopic disorder, R s 3> r LOl can be achieved, 
for example, in experiments with 4 He films absorbed on a 
Mylar substrate when the torsional oscillator technique is 
implementedi 23 i 26 In these cases, typical u> ~ 27rx 1500Hz 
leads to r^ ~ 5 /im. Mylar contains bubbles and flakes of 
sizes 1— 100/xm and average separation of 10— 150^m, i.e., 
R s ~ 100 /im. Macroscopic variations of the substrate 
potential result in the inhomogeneous broadening of the 
KT type superfiuid transition. Similar to our suggestion, 
frequency dependence of the inhomogeneous broadening 
of the jump of the superfiuid density and the dissipa- 
tion peak was successfully studied by a two-torsional- 
oscillators technique^ 

The fact that in our simulations we relate rough sur- 
face of the CaF 2 substrate of Ref. Q] with the spatially 
varying potential implies that helium film thickness is 
modulated by the microscopic surface roughness. This 
is true for the case of the superfiuid transition at high 



IV. RELATED PHYSICAL SYSTEMS 

Our research is centered around the KT-type 2D finite 
temperature superfiuid transitions in 4 He films. Phase 
transitions in thin-film superconductors^ Josephson 
junction arrays^ weakly interacting 2D Bose gaseS ) 21 i 31 
2D gas of spin-polorized atomic hydrogen on liquid- 
helium surface j 3 - 2 . and Bose-Enstein condensates loaded 
on a 2D optical lattice 3 - 3 - belong to the same universal- 
ity class. Our analysis of the role of diagonal disorder in 
4 He films can be qualitatively extended to some of these 
cases. 

Microscopic details are irrelevant for the study of uni- 
versal properties of the second order KT phase tran- 
sition characterized by the diverging correlation length 
scales ] 13 ! 16 For example, universal jump |T]) of the areal 
superfiuid density is innate for both discrete lattice type 
and continuous geometry of experimental realizations of 
2D bosonic systems with diagonal disorder and controlled 
value of the chemical potential. Instead of 4 He film ab- 
sorbed to rough surface with the chemical potential con- 
trolled by the pressure in the sample cell^ we could have 
considered as well an atomic gas in an optical lattice with 
a given number of atoms per site and a varying trapping 
potential, or a Josephson junction array (or a supercon- 
ducting granular film) with an external voltage applied 
to the ground plane and random gate voltages induced 
by trapped charge in the substrate. 34 Critical behavior 
in all mentioned cases falls within the framework of the 
lattice bosonic Hubbard model ([2]). The same finite-size 
analysis as in Sec. Ill Bl can be used to resolve the univer- 
sal KT jump However, microscopic details determine 
how strong will be finite-size effects. Thus, similar to our 
study targeting 4 He films absorbed to rough surfaces in 
Ref. [l|, separate simulations are required for the particu- 
lar experimental situation to reveal (or disproof) primary 
importance of finite-size effects as a source of experimen- 
tal smearing of the discontinuous KT transition. 

Along with diagonal disorder, off-diagonal disorder 
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can take place in 2D systems that exhibit finite tem- 
perature KT transition. Off-diagonal disorder can be 
present in the form of site or bond dilution disorder in 
Josephson junction arrays 3 -^, or as nonmagnetic disor- 
der in planar symmetry spin models 3 ^, or random gauge 
disorder— (such as disorder in the phase of the complex- 
number tunnelling amplitudes for bosonic systems), etc. 
Our study of non-universal corrections to the finite-size 
broadening of the KT transition produced by diagonal 
disorder cannot be directly related to the case of off- 
diagonal disorder. 2D XY models are used for descrip- 
tion of aforementioned disordered cases. But an XY 
model simply does not admit conventional diagonal dis- 
order. Formally within the bosonic Hubbard hamilto- 
nian (JSJ) the XY limit is obtained by writing the last two 
terms as U/2 [hi — {p, — Vi)/U) 2 and setting simultane- 
ously U — * oo, and /i — > oo. In this limit average number 
of particles per site is fixed, n = fi/U, and there are no 
local density fluctuations at all, i.e., any finite diagonal 
disorder vi in /j, is ignored with respect to the XY models. 

On a passing note, finite-size effects and disorder will 
play an important role in the observation of the supercon- 
ducting transition in 2D (the KT scenario is debated now 
for low-T c granular films, 38 high-T c films^ and Joseph- 
son junction arrays42.) However, we do not continue our 
discussion of both neutral superfluid and superconduct- 
ing [charged superfluid] systems due to differences in the 
corresponding measurement techniques, and due to addi- 
tional magnetic field screening specific for superconduct- 
ing systems.— 

V. DISCUSSION AND CONCLUSIONS 

We have studied a 2D bosonic system with diagonal 
disorder in the vicinity of the finite temperature super- 
fluid transition. We found that the disorder dramatically 
increases finite-size smearing of the universal jump of the 
superfluid density ((T|). And, an analysis of its finite-size 
scaling is required to resolve faithfully the underlying KT 
scenario of the superfluid transition. 

Our simulations consistently explained difficulties of 
identification of the KT transition by QCM frequency 
shift in the experiments of Ref. [l| with 4 He absorbed to 
microscopically rough CaF2 surfaces. In those experi- 
ments QCM measurements were carried out at a fixed 
frequency w. Thus, for a given disorder strength, the 
effective finite length scale set by the vortex diffusivity 
length r w oc l/y/u> was fixed. Without finite-size scal- 
ing it was not possible to extract the universal jump of 
the superfluid density (fTJ) in the case of strongly disor- 



dered substrates. We speculate that operation of QCM 
crystals at different harmonics of their fundamental res- 
onance will be required. Note that QCM measurements 
of the superfluid transition of 4 He films on a flat sub- 
strates at two different harmonics of fundamental reso- 
nance frequency were carried out. Finite-size broadening 
that increased at higher frequencies was observed^ An- 
other study of the u dependance of the KT transition in 
2D 4 He films used an ultrasonic technique^ 3 - For a given 
sample the superfluid transition was relatively broad at 
high frequencies with the universal jump of the superfluid 
density smeared toward the normal phase, but sharpened 
as the frequency was decreased (correspondingly, the po- 
sition of the dissipation peak was shifting). Our simu- 
lations (Sec. IIII B[) predict the same scenario for a given 
disorder strength as r u is increased. 

The experiments of Ref . [l] were our original motivation. 
Due to the noise in the measured data we sought only a 
qualitative description of those experiments. Thus, we 
used a lattice Hubbard model with on-site potential dis- 
order to describe helium films and discrete space and 
imaginary time worm algorithm to carry out classical 
Monte Carlo simulations (Sec. [TT|- Universal properties 
of the superfluid transition were preserved, but one-to- 
one matching of the experiments and the simulations can- 
not be achieved in this kind of simulations setup. Exact 
ab initio numerical simulations of 4 He film with realistic 
substrate potentials and geometry is possible with exten- 
sion of the used worm algorithm 17 for continuous-space 
path Monte Carlo simulations^ Another aspect is that, 
with regard to the experiments of Ref. ,l| on 4 He films, we 
have considered the case of a 2D bosonic system driven 
to the point of the superfluid transition by the increase 
of boson density (i.e., chemical potential was the con- 
trol parameter). We expect qualitatively similar results 
for broadening of the KT transition if, instead, at a con- 
stant boson density the system will be brought across the 
superfluid-normal phase boundary by an increase of the 
temperature (see Fig.QJ. 
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